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We present color-magnitude diagrams of globular clusters for models with self-enrichment and 
pre-enrichment. The models with self-enrichment turn out to have two or more main sequence 
turnoff points in the color-magnitude diagram if the fraction of mass lost by the globular cluster 
under supernova explosions docs not exceed 95-97%. The models with pre-enrichment can have 
£^ . only one main sequence turnoff point. We argue that the cluster w Ccn evolved according to a 

| self-enrichment scenario. 

1 INTRODUCTION 

Globular clusters (GCs) belong to the oldest stellar population of the Galaxy. The ages of the 
oldest GCs (13-14 Gyr) suggest that they have formed as distinct structures at the very beginning 
of the Galaxy formation, possibly preceding formation of the field stellar population. Therefore, it 
is natural to assume that heavy elements in GCs have been synthesized during the process of stellar 
evolution in the GCs themselves, and were not brought in with the gaseous clouds from which the 
clusters formed. This point of view is discussed recently as the self-enrichment hypothesis, contrary 
to the pre-enrichment hypothesis, which assumes that the abundances of metals currently observed in 
GCs are essentially the same as their values in the gaseous clouds of proto-GCs; in other words that 
the metals have been produced by an external stellar population [1]. Note that the pre-enrichment 
scenario seems to be required for at least some GCs (such as the young group of GCs). This follows 
from the bimodality of GCs, i.e. from the fact that they form two groups with different metallicities, 
as was first noted in [2] (for later discussions see [3]). We stress that the concept of self-enrichment 
does not exclude a possibility that proto-GCs gaseous clouds have already contained some metals 
produced by the very first Population III stars, though it suggests that their amount was negligible 
compared to the observed values. 

Neither of these hypotheses can currently be rejected or accepted confidently based on observations. 
Arguments exploiting expected mass-metallicity correlation [4] cannot be applied to GCs directly due 
to their low binding energy [5]. Furthermore, some evolutionary scenarios for GCs predict an anticor- 
relaton between their mass and metallicity [HE]: Z oc M _1 . A weak anticorrelation between mass and 
metallicity discovered in [5] for old GCs might argue in favor of the self-enrichment hypothesis, but is 
not entirely convincing due to the small sample studied {y~N — 4.5); in addition, this anticorrelation 
is unstable to effects related to the decrease in the masses of GCs due to evaporation [7], gravitational 
encounters [8], and tidal forces [9]. Under these circumstances, it seems necessary to search for possible 
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manifestations of self-enrichment in the color characteristics of some particular GCs, and for possible 
expected differences between the characteristics of stars in old and young GCs. 

In this paper, we present the results of modeling of the fine structure of the color-magnitude 
diagrams of GCs in the vicinity of their main sequence turnoffs, and describe the dependence of 
this structure on enrichment scenarios. Recently published detailed observational studies of the GCs 
uj Cen |10] and NGC 2808 [UJ provide information about the star formation history (SFH) and 
enrichment of these GCs. A comparison of our modeling results with the observed color-magnitude 
diagrams for uj Cen and NGC 2808 provides support for the self-enrichment scenario for both uj Cen 
and NGC 2808. 

uj Cen is distinguished among GCs in its unusually high mass (M ~ 3 • 1O 6 M ) and the large 
scatter of stellar metallicities |12j . The latter seems to provide evidence in favor of self-enrichment 
in this cluster [13]. Furthermore, as the binding energy of uj Cen does not substantially exceed the 
binding energies of most GCs in the Galaxy p3], this suggests that at least some fraction of the 
Galactic GCs were enriched in metals by nucleosynthesis occurring inside the clusters themselves. 
The absence of stars with metallicities [Z] < — 2 in uj Cen does not necessarily require pre-enrichment 
to this metallicity as noted in [12j . but may provide evidence for triggering of star formation in uj Cen 
by activity of supernovae in the cluster (see [6] for discussion). 



2 DESCRIPTION OF THE MODEL 

The parameters of the GCs were computed using a single-zone model [6] based on the standard 
system of equations describing evolution of the gas mass and the elemental abundances, and a single- 
zone imitation of dynamics associated with the energy released by supernovae [151 IT6| \T7\ I18j. 

2.1 Chemistry 

The first model equation is the conservation of the gas mass: 

_^£ = -^(t)+ J ij(t-Tj^(M,t-Tj(M-M r )dM-M° ut + Mi n , (1) 

where M g is the mass of gas in the system, M r (M), the mass of stellar remnants, tp(t), the star 
formation rate (SFR), <p(M,t), the initial mass function (IMF) at time t, M m i n and M max are the 
minimum and maximum masses of new-formed stars, r u is the lifetime of a star with mass M, and 
Mout an( j j^in are ra ^ es Q f ejection and accretion of gas by the system, respectively. 
We used the standard Schmidt law pH [T?J for the SFR: 



ip(t) = fp n V, (2) 

where / is the star formation efficiency (SFE), p, the average density of the gas, and V, the volume of 
the stellar system. Schmidt [20] adopted the value n = 2. The quadratic dependence follows from star 
formation models with self-regulation, in particular, the model in which star formation is regulated 
by the ionization of gas by UV photons from massive stars [21]. The efficiency / is determined 
by the specific mechanisms providing a power-law dependence of the SFR, ip ~ p n , and in general 
determination of / is a separate problem. A good estimate for the SFE in the Schmidt law is probably 
provided by / = 2 • 10 7 sm 3 g~ 1 s~ 1 , as adopted in [191 IT?] . 
We used a Salpeter IMF [22]: 

ip (M ) ~ M" 2 - 35 (3) 

with the masses in the range from M m j„ to M max . 

The equation describing time variations of the radius of the gaseous component R g is 
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dR g _ Rlje R SN (t) R g 



(4) 



dt 3GM GC M g R g 6r d ' 

where M GC is the total mass of the GC, R* is the radius of the stellar component, while R g is the 
gas component radius, G, the gravitational constant, 7, the efficiency of the transformation of the 
supernova explosion energy into the energy of turbulent gas motions, £0, the explosion energy, r^, 
the timescale for energy dissipation in cloud collisions, and R SN (t), the rate of type II supernovae. 
Equation ([4]) is a spherically symmetric version of the dynamical equation written fist by [16] for a 
flat geometry. 

The type II supernovae rate R SN (t) is defined by the integral 

R SN (t)= J tp (t — t m ) (p (M, t — t m ) dM, (5) 

Mi ow 

where Mi ow = 8M Q is the minimum mass of a presupernova. The evolution of the chemical composition 
of a GC is described by the equation of mass conservation for a particular element: 
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where Zi(t) is the relative abundance of element i at time t; Xi(M, t — t u ) and X- a (M, t — t m ) are the 
masses of element i ejected by evolved stars and by binaries of mass M that were born at time t — r M ; 

and M™. describe the exchange of heavy elements with the intragalactic medium through wind 
and accretion, respectively; M\°™ = 1.5M and = 8M Q are the minimum and maximum masses 
of stars in binaries; (5 is the fraction of binaries in the mass interval 2M l j™-2M^ ; and /(//) is the 
mass distribution of stars in binaries (where \x is the ratio of the mass of the lighter component and 
the mass of the binary, which varies from fJ, m in = max(M#7M, 1 - M^/M) to 0.5). 

The first term in (0) describes the decrease of Z{ in the interstellar medium due to star formation; 
the second term corresponds to enrichment of the interstellar medium with ith element by stars with 
masses in the range ^M m j n , 2Mj° w ^ ; the third term describes the contribution from single stars with 

masses M = (2Mj", 2M"^); the fourth term describes the input from Type la supernovae in binaries 
with masses ranging from 2M{° tu to 2M^; the fifth term is the contribution to Z{ due to stars more 
massive than 2Mf^; the last two terms describe the exchange of heavy elements between the cluster 
and the external medium. 
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In the present study, the value of (M, t — t m ) for stars with masses M < Mi ow is determined 
by equation 

Xt (M, t-T M ) = Z i (t- r M ) (M - M r (M)) . (7) 

The yields of heavy elements by Type II supernovae we computed with using the results of Woosley 
and Weaver [23], while for Type la supernovae we used the yields from [23]. The latter yields are 
valid for solar metallicity only, but the yields of Type la supernovae do not depend substantially on 
the metallicity; the mass distribution of stars in binaries we took in the form [25]: 

/M = 2V (8) 

2.2 Photometry 

At time t the system retains the following number of stars with the masses from M to M + dM born 
at time from tf to t' + dt': 

dN = ^(t')ip(M, t')Q(r M - r)dMdt', (9) 

where @(x) is the Heaviside function. For a given stellar mass and age r = t — t', the stellar luminosity 
L and effective temperature T can be found from a database of evolutionary tracks. For this purpose 
we used the Padova data base [26] [271 EH ESI ED] supplemented with data for low mass stars from [3T] . 
For the color-magnitude diagram we have utilized the widely used library of stellar spectra described 
in [32[ [33] with the data for stars of effective temperatures T < 50000-fT, and the supplementary 
library [34] with the data for hot stars T > 50000^. If one defines the spectrum of a star as 

F(X;L,T,g)=L-f(X;T,g), (10) 

where L is the luminosity in solar units, f(X;T,g), the spectrum normalized to the solar spectrum, 
and g, the surface gravity, the absolute magnitude is given by 

oo 

m = -2.5 lg/ F{\)$(\)d\ + m cal]h 
o 

oo 

= -2.51gL-2.5lgJ f(\)<S>(\)d\ + m calih , (11) 
o 

where m ca \fo is the calibration constant of the filter corresponding to the given waveband, with the 
filter sensitivity curve <&(A) [35]. 

3 RESULTS OF MODELING 

Earlier [6] we have considered and theoretically motivated a model of chemical evolution of GCs 
within the hypothesis of self-enrichment with a varying IMF: the initial stage of star formation was 
characterized by a "top-heavy" IMF shifted to M > 8Mq, while the transition to a second stage of 
star formation with a "normal" IMF occurred after heavy elements started to dominate in radiative 
cooling of the proto-cluster gas. Note, that the possibility for the first stars forming of pristine matter 
to have predominantly high masses is confirmed recently in numerical experiments (see, [36]). 

Here, we present the results of computations of color-magnitude diagrams V-(B-V) for a similar 
model. Figure 1 shows the V-(B-V) plane for a GC model with the mass of M = 1O 6 M0 that 
undergoes two episodes of star formation: the first with a Salpeter-like IMF, but shifted to higher 
masses, so that the minimum mass of the stars formed is M m j n = 8M Q ; the maximum mass of stars 
is assumed M max = 100M Q for all cases. The SFE in the first stage is taken f\ = 2.6 • 10 5 cm 3 g _1 s _1 , 
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and the transition time to the second stage (with a normal IMF) is t = 10 Myr. We assumed for 
the second episode of star formation M m i n = O.IMq and /2 = 2 • 10 7 cm 3 g _1 s _1 . In order to avoid 
overproduction of heavy elements in supernovae explosions, we assumed that 95% of the matter ejected 
by supernovae is swept up after i = 10 Myr. In these assumptions two main sequence turnoffs appear 
clearly in the color-magnitude diagram. 

With the model parameters shown above, the presence of several turnoffs is a natural consequence 
of the evolutionary processes of the stellar system, if transformation of the IMF occurs as assumed in 
[6] when the metallicity exceeds some critical value: during the transition to the standard IMF the 
SFR attains the peak value ip ~ 0.6 M & /yr due to increase in the SFE by two orders of magnitude. 
The SFR remains near its maximum level for 1 Myr. During this time, a significant fraction of stars 
(~ 30%) form from gas that has already been enriched by the first supernovae. The stars have formed 
in this period have metallicities [Z] ~ —2 and concentrate mainly toward the left side of the diagram. 
A fraction of them fall by the present into the region of the left subgiant branch, as seen in Fig. [TJ After 
the transition to the "standard" IMF, the metallicity of the stars in the system continues to increase. 
At time ~ 2 • 10 8 yr, the metallicity attains a practically constant asymptotic value [Z] ~ —0.8. In 
further continuing star formation, an additionanl population of stars is born in the cluster, now with 
a higher metallicity, close to the asymptotic value (~ 10%). 

By varying the free parameters of the model, such as the time for transition to the standard 
IMF, star formation efficiencies in the first (/i) and second (/2) episodes, the time when the matter 
from supernovae explosions begins to leave the system, and the fraction of lost matter, one can 
obtain even more complex color-magnitude diagrams in this scenario. An example is shown in Fig. [2 
which corresponds to SFE on the second stage an order of magnitude higher than before (/2 = 
2 • 10 8 sm 3 g~ 1 s~ 1 ), with all other parameters fixed. We stress that the transition from the diagram 
with two turnoff points to the one with three turnoff points is reached in this case by increasing 
only one parameter ■ - j% — without necessty to regulate "by hands" the sequence of several star 
formation bursts. In turn, the SFE may be determined by a number of factors, which at present 
cannot be estimated confidently. For instance, in a model in which the star formation is self-regulated 
by ionizing stellar radiation [21j . the coefficient /2 may be critically determined by inhomogeneity of 
the interstellar gas after the first burst of star formation and the fraction of ionizing radiation able to 
escape the proto-cluster. 

An increase in the SFE results in a higher than in the latter case number of stars in the left side 
of the diagram. On the other hand, it leads to a faster depletion of gas, and therefore to a higher 
enrichment in metals {[Z] ~ —0.3). A fraction of the stellar population born during the period when 
the metallicity of the system is maximum reaches currently the extreme right of the subgiant branch. 
Later on, the gas metallicity decreases due to mass loss of evolved low- and intermediate-mass stars, 
and attains a constant level, corresponding to the one of stars currently located in the middle subgiant 
branch. 

In the case considered here the presence of several turnoff points is a consequence of natu- 
ral evolution of the stellar system, and does not require stronger assumptions about the nature 
of the star formation, or the existence of a stellar population with unusually high helium abun- 
dance (Y ~ 0.4 [371 ESI EHl HQ]), which would require an increase in the metal abundance by 
AZ ~ AY/ 5 ~ 0.08. The diagrams in Figs. Q] and [2] are plotted with accounting the observa- 
tional errors. We used an approximation of the data from [H] to calculate the magnitude dependence 
of the rms errors. 

The presence of several turnoff points indicates also that the main sequence itself has a complex 
structure. However, it is difficult to distinguish this structure from such kind of color-magnitude 
diagrams, mainly because of large observational errors for low-mass stars. In order to analyze the 
main sequence, we calculated several color distributions of stars along the "cuts" corresponding to a 
fixed value of V as suggested in [10] . Figures [3] and [4] show the dependences of the relative numbers of 
stars for three values of V. Figure [3] corresponds to the color-magnitude diagram shown in Fig. Q] and 
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Fig. U] to the diagram shown in Fig. [5J In these figures left panels show the distributions accounting 
the Gaussian smearing described above; for comparison, we show in the right panels the distributions 
obtained for the same model without contribution from observational errors. These figures show that 
even though the system does possess several main sequences, the distinction between them will not 
always be revealed in the observational color-magnitude diagrams. 

A group of stars located on the continuing of the main sequence to V ~ 1, in the region occupied 
by so-called blue stragglers, is clearly seen. In our models, these are the youngest stars, born in late 
stages of evolution of the GC. For instance, the group of stars in Fig. [2] that can be identified with 
blue stragglers is characterized by lower metallicities ([Z] ~ —1) than the stars located in the right 
side of the diagram ([Z] ~ —0.3). This relatively small group of blue stragglers is born in the model 
shown in Fig. [2] at late evolutionary stages of stellar system when the metallicity of the gas begins to 
decline due to mass loss by low-mass stars with lower metallicities. Since the gas mass remained in 
the system by this time is small, the mass income with lower metal abundances becomes sufficient to 
significantly decrease the gas metallicity [6]. This identification is not acceptable in reality, however, 
since blue stragglers typically have anomalous chemical composition, which probably indicate their 
connection to active mass exchange in binaries. Note that there are no stars in the region of the blue 
stragglers in the color-magnitude diagrams if the model assumes that star formation does not occur 
at ages exceeding 1 Gyr. 

3.1 A Model with Pre-Enrichment 

Models with self-enrichment assume a longer evolution of the system than models with pre-enrichment. 
Indeed, models with pre-enrichment suggest that the bulk of metals has been already injected into the 
gas by some external sources on the proto-cluster stage. Therefore, over the entire subsequent cluster 
evolution only a slight enrichment is allowed, which does not significantly exceed the amount of the 
pre-enrichment. In turn, this means that in models with self-enrichment star formation has to occur 
in a mode that prevents the birth of stars from newly enriched matter. In other words, in such models, 
newly enriched matter must be swept out of the cluster before it will be able to cool radiatively and 
give rise to new generations of stars. A sufficient condition for such a mode to occur is a rapid star 
formation on early stages, which results after several Myr in a large energy release enough to prevent 
the enriched gas from cooling. A rough estimate of this condition can be obtained by assuming that 
the supernovae energy rate exceeds the total rate of radiation losses in the gas: 

u sn E sn > A(T)n 2 V, (12) 

where u sn and E sn = 10 51 erg are the supernova rate in a GC, and the supernova explosion energy, 
A(T) ~ 3 • 10 -23 erg sm 3 s _1 is the rate of radiative losses at gas temperatures behind the shock 
T ~ (1 — 3) • 10 5 K, n is the number density of the gas and V is the volume of a GC. For n ~ 10 3 ~ 3 
and R ~ 20 pc, the estimate (fT2j) gives v sn ~ 3 supernovae per year, which is equivalent to a SFE 
/ ~ 2 • 10 8 sm 3 g~ 1 s~ 1 an order of magnitude higher than commonly assumed value. This is obviously 
an overestimate, since the binding energy of globular clusters is low, and a substantially lower energy 
release is sufficient to remove the gas from the cluster. Indeed, it is easy to show that, even for the 
SFE assumed above, / ~ 2 • 10 7 sm 3 g~ 1 s~ 1 , the number of supernovae that would explode in the 
cluster over time sufficient for the volume filling factor of the hot gas to equal one is N sn ~ 600, whose 
total energy release exceeds the binding energy of a cluster with mass M = 1O 6 M0 and initial radius 
R = 20 pc by two orders of magnitude. Below we will adopt precisely this SFR for the model with 
pre-enrichment, / = 2 • 10 7 sm 3 g~ 1 s~ 1 . 

Figure [5] shows the results of computing the color-magnitude diagram for a GC with a mass of 
M = 1O 6 M which is pre-enriched to metallicity [Z] = —1.5 and begins its evolution with radius 
R = 20 pc and the SFE given above. Similar to the previous figures, Fig. [5] is plotted with accounting 
a Gaussian smearing. An important feature of this model is that, when the first supernovae begin to 
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eject heavy elements at about t > 10 7 yr, the bulk of the matter with the initial chemical abundances 
is already contained in stars. Thus if we assume that under these circumstances the gas ejected by 
supernovae is easily removed from the system, the number of stars with metallicities exceeding the 
initial value will be negligibly small. Figure [5] shows the results for a model GC in which practically 
all stars have the initial metallicity. 

It is obvious that a similar diagram can also be obtained in a model with self-enrichment if the 
parameters of star formation are such that the separation (in B — V) between the main sequences is 
smaller than the observational errors. In this case, it is practically impossible to distinguish between 
the two populations on observational color-magnitude diagrams. As an illustration of this situation, 
Figures [6] and \7\ present a color-magnitude diagram and color distributions of relative numbers of 
stars, respectively. In this model, the first burst of star formation has M m i n = 8M , M max = IOOMq, 
and fi = 2.6 • 10 5 sm 3 g _1 s _1 ; the transition to the second burst with a "standard" IMF occurs 
at t = 10 Myr; the second burst of star formation has M m i n = 8M & , M max = 100 M & , and fa = 
2 • 10 7 sm 3 g~ 1 s~ 1 , with 99.7% of the matter ejected by supernovae swept out of the cluster after 
t = 10 Myr. In this case, the absence of multiple turnoff points is a consequence of a larger fraction 
mass loss than in the model shown in Fig. [U Note that the model computation produces two stellar 
populations which correspond to low metallicities [Z] = —2 and -1.6; this makes it even more difficult 
to distinguish between the populations, due to the proximity of the isochrones for low metallicities. 

Thus, the absence of multiple turnoff points in observational color-magnitude diagram does not 
necessarily mean that matter of the corresponding GC was pre-enriched in heavy elements, since the 
separation that arises in models with self-enrichment may be smaller than the errors in the colors. In 
self-enrichment scenarios, it may be easier to distinguish stars with different metallicities within more 
massive GCs. In fact, we would expect a decrease of the metallicity of the main fraction of stars with 
increasing cluster mass [4j [6] . At the same time the later generation of stars has as a rule metallicities 
close to the asymptotic value \Z\ ~ —0.8, which only weakly depends on the cluster mass. Observing 
this separation for GCs with masses of M ~ 10 5 Mq requires an increased accuracy for magnitudes 
of Am ~ 0.03 — 0.07. The proposed interpretation of the color-magnitude diagram of co Cen as 
having multiple main sequence turnoff points with several distinct stellar populations is justified by 
spectroscopic measurements of metallicities in co Cen. Several main sequences are also observed in 
NGC 2808 [43]. With accounting the fact that spectroscopic metallicities of different branches of the 
main sequence in NGC 2808 are nearly equal, Piotto et al. [32] suggest a scenario where the multiplicity 
of the main sequence is connected with different helium abundances (Y ~ 0.24 — 0.4). Thus, the two 
cases — co Cen and NGC 2808 — seem to represent two possible scenarios for the evolution of GCs: 
one produces several populations with different metallicities, while another produces populations with 
similar metallicities but very different helium abundances. From this point of view in order to clarify 
the origin of this multiplicity, detection of multiple main sequence turnoff points in other Galactic 
GCs and a precise spectroscopic determinations of the iron abundance will play crucial role. 

4 CONCLUSIONS 

We have applied a simple single-zone model to analyze color-magnitude diagrams for globular 
clusters. Our results show the following. 

1. The presence of two or more main sequence turnoff points in the color-magnitude diagrams of 
GCs can be naturally explained in models with self-enrichment, with an initial stellar mass function 
that is initially shifted toward high masses, but then changes to a normal IMF. 

2. This represents an argument in favor of the self-enrichment of GCs with multiple turnoffs (such 
as co Cen) in metals due to internal processes within the GC. 

3. On the other hand, the uncertainties in the luminosities and magnitudes of stars in GCs are still 
fairly large, and do not always enable resolution of the main sequences in color-magnitude diagrams 
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into subgroups. In other words, the absence of multiple turnoff points does not necessarily require an 
external (pre-) enrichment in metals. 
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Figure 1: Color-magnitude diagram for bimodal star formation. The initial star formation episode 
has M m i n = 8 M & and f± = 2.6 x 10 5 sm 3 g— Is -1 , and lasts for 10 7 yr; the second star formation 
episode has M m , n = 0.1 M and fi = 2 x 10 7 sm 3 g _1 s _1 . At t = 10 7 yr, 95% of the matter ejected 
by supernovae has left the cluster. 
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Figure 2: Same as in Fig. 1 for fi = 2 x 10 8 sm 3 g 1 s . 
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Figure 3: Number of stars as a function of the color corresponding to cuts of the color-magnitude 
diagram shown in Fig. 1 for V = 5, 6, 7. Diagrams (a), (b), (c) correspond to the plotted histogram, 
taking into account the observational errors, while diagrams (d), (e), (f) correspond to the plotted 
histogram without accounting observational errors. 
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Figure 4: Same as in Fig. 3 for analogous cuts of the color-magnitude diagram presented in Fig. 2. 
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Figure 5: Color-magnitude diagram for the model with primordial self-enrichment [Z\ n = —1.5; the 
SFR is specified to be zero at t = 10 s yr. 
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Figure 6: Same as in Fig. 1, but with the fraction of matter ejected by supernovae and lost by the 
system at time t = 10 7 yr being 99.7%. 



12 



0.2 
0.1 
0,0 
0.6 
0,5 
S 0,4 
? 0,3 
0,2 
0.1 
0.0 



V=4 



V=4 



.4 



d) 



V=5 



b) 



V=5 



e) 



V=6 



V=6 



^^^^^^^^ 



0,0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 0,0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0 



B-V 



Figure 7: Same as in Fig. 3 for the analogous cuts of the color-magnitude diagram presented in Fig. 7. 
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